Antikaons and higher order couphngs in relativistic-mean field study of neutron stars 



Ncha Gupta and P. Aruniugam 
Department of Physics, Indian Institute of Technology Roorkee, Roorkee 



247 667, India 



m 
o 

(N 



> 
00 

ON 

o 

m 



X 



We investigate the role of higher order couphngs, along with the condensation of antikaons {K~ 
and on the properties of neutron star (NS). We employ extended versions of the relativistic 

mean-field model, in which kaon-nucleon and nucleon-nucleon interactions are taken on the same 
footing. We find that the onset of condensation of K~ and K'^ highly depends not only on the 
strength of optical potential but also on the new couplings. The presence of antikaons leads to a 
softer equation of state and makes the neutron star core symmetric and lepton-deficient. We show 
that these effects strongly influence the mass-radius relation as well as the composition of neutron 
star. We also show that the recently observed 1.97±.04 solar mass NS can be explained in three 
ways; (i) a stiffer EoS with both antikaons, (ii) a relatively soft EoS with K~ and (iii) a softer EoS 
without antikaons. 



PACS numbers: 26.60.-c,26.60.Kp, 13.75. Jz,97.60.Jd 



I. INTRODUCTION 

Neutron stars (NSs) are fascinating objects, attracting 
strong appeal as probes into the understanding of many 
areas of physics. The interior of a NS, where the density 
is very high, provides opportunities to apply and test the 
concepts of nuclear physics in order to elucidate proper- 
ties such as NS mass and radius. These properties depend 
on the microscopic nature of the matter at high densities 
and can be corroborated with astrophysical observations 

In the recent years, several new observations have pro- 
pelled several theories anda variety of models in a bid 
to understand the properties of NS frl-fisj . In NS mat- 
ter, while considering the conservation of only charge and 
baryon numbers, antikaons, hyperons and quarks can ap- 
pear inside the NS by a strangeness changing process. 
At intermediate densities, pions and antikaons {K~ and 
K^) are most likely to condense. In a vacuum, pions are 
lighter than kaons, but this situation may be reversed 
in the dense medium due to strong interactions between 
mesons and nucleons [l^ [l^ • Hence antikaon condensa- 
tion becomes more important in the intermediate density 
range and affects the mass and radius of the NS signifi- 
cantly. 

As we approach the interior of the NS, density in- 
creases. The excitation energy of antikaons (strangeness 
= — 1) decreases with density and hence at a sufficiently 
high density the antikaons are favoured to condense. 
Many theoretical descriptions have considered only the 
condensation |18l - [24{ which is supposed to be domi- 
nant as it occurs at a relatively low density. With the on- 
set oi K~ condensation, n p-\-K~ is the most preferred 
process, and hence the proton fraction rises dramatically 
and even exceeds the neutron fraction at higher densities 
. With the onset of K° condensation this scenario will 
change completely. There will be a competition between 
the processes N N + and n p + K~ , resulting in 
a perfectly symmetric matter of nucleons and antikaons 
inside the neutron stars [l9| . 

A detailed study of the presence of antikaons, initially 



proposed by Glendenning and Schaffner-Bielich [I^ , with 
relativistic mean-field (RMF) models has been carried 
out in [25I - I27I . In most of these previous works, the 
RMF paramctrizations used were not tested rigorously 
in the case of finite nuclei. The recent paramctrizations 
discussed in this work have been carefully developed and 
tested over the nuclear chart by explaining several nu- 
clear properties [28^30} and the experimentally deter- 
mined EoS for symmetric nuclear matter [lo| . Hence, we 
are using more reliable models for the EoS which turn out 
to be softer (around the saturation density) than those 
used in the previous works. It is quite well-known that 
a softer EoS could lead to a lesser contribution from an- 
tikaons [2^ . We have shown in our earlier work that 
the EoS from recent parameterizations, with higher order 
couplings, is not too soft to neglect the role of antikaons. 
In such a case we have shown that the mixed phase (of 
kaonic and nonkaonic phases) will not appear due to the 
softer EoS. It has to be noted that it is not only the stiff- 
ness of EoS but also that of the symmetry energy that 
is crucial in determining the onset of antikaon conden- 
sation; hence their role in modifying the properties of 
NS. Some of the recent parameterizations (e.g. FSUG- 
old) yield very soft symmetry energies and its effect in 
altering onset oi and hence the properties of NS are 
discussed in our earlier work . This study introspects 
our earlier conclusions with the inclusion of K'^ which is 
a more realistic case [loj . 

For our calculations, we consider the effective field 
theory-motivated relativistic mean-field model (E-RMF) 
jsif . In this model, the idea of renormalizability is aban- 
doned and the effective Lagrangian is expanded in powers 
of fields and its derivatives at a given order, with all the 
non-renormalizable couplings consistent with the under- 
lying symmetries of QCD. In short, one can say that the 
E-RMF model comprises the standard RMF plus a few 
additional couplings. The RMF terms dominate at low 
density while the additional couplings dominate at high 
density. Without forcing any change in the parameters 
initially determined from a few magic nuclei [3l|, the E- 
RMF calculations explain finite nuclei and infinite matter 
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in a unified way with a commendable level of accuracy in 
both cases [l^. It is interesting to see how the E-RMF 
description of NS changes with the inclusion of antikaons, 
which is the central subject of this work. 

In section |lll we describe the Lagrangian, field equa- 
tions and the expression for energy density followed by a 
discussion on the constraints for the antikaon condensa- 
tion, with the parameters used. Our results and discus- 
sions are presented in section Hill which is followed by the 
conclusions drawn from present work. 



II. THEORETICAL FRAMEWORK 

The effective Lagrangian for the extended RMF mod- 
els, after curtailing terms irrelevant to nuclear matter, 
can be written as 
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For the FSU2.1 model, the above Lagrangian has an addi- 
tional term K^glRf^.R^'glV^V^' . The symbols g^, g^i, gp, 
'*3j rji, 772, rip, A-u and Co denote the various coupling 
constants, cr, V^, and ijj denote the scalar, vector 
and isovector meson fields and the nucleon field, respec- 
tively. TOcr, "rriuj, and nip are the corresponding meson 
masses and m„ is the nucleon mass. More details of the 
Lagrangian are explained explicitly in [ij, |3l| . 
The Lagrangian for the antikaon part reads 



Ck = D*K*D^'K 



'-K*K, 



(2) 



with K = or K° and is added to the E-RMF La- 
grangian. The scalar and vector fields are coupled to 
antikaons in a way analogous to the minimal coupling 
scheme [l^ via the relations 

= rriK — gaK<^ and (3) 
Dp = dp + igujKVp + igpKT3.Rp, (4) 



where niK stands for the antikaon's mass {mj^- = 
rrij^o = 495 MeV). Note that in the mean- field approx- 
imation, the fields cr, Vp and Rp are replaced by their 
expectation values cr, Vq, and i?o, respectively. In the 
presence of antikaons the coupling constants correspond- 
ing to these fields arc represented by g^K , g^K and gpK ■ 
Energy relations for the antikaons {K~ ,K^) arc 



^K-,K'> — "^K — gaK<^ — guiKVo =F gpRRo, 



(5) 



where =F sign represents the isospin projection of an- 
tikaons K~ and K'^ respectively. The expression for the 
energy of antikaons is linear in the meson field and repre- 
sents that, with the increase of density, the energy of an- 
tikaons will decrease. The above expression also suggests 
that antikaon condensation is significantly infiuenced by 
the rho meson field or vice-versa. 

In the presence of antikaons, the meson fields are given 

by 
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rnpRo = x3p(Pp - Pn) ~ lip mpRo 
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-gpK{PR- - Pro), 
where ps is the scalar density given by 
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and 7 is the spin-isospin degeneracy factor and is equal 
to 2 (for spin up and spin down). 

The densities of antikaons can be written as, 

Pr-,ro = '2'{^R-,R0 + gujRVo±gpRRo)K*K. (10) 

In NS matter only baryon number and the charge number 
arc conserved, hence the constraints involving chemical 
potentials and baryon densities can be written as 
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The energy density can be written as 
e = En + £r-,ro, 



(11) 



(12) 



where £nJs the energy density of nucleon phase as given 
in Ref. [24| . The energy density contributed by antikaon 
condensation is 



,R0 = m^iPR- + Pro)- 



(13) 



Unlike the energy density, pressure is not directly affected 
by the inclusion of antikaons in an s-wave condensation 
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[igj . The conditions for onset of antikaons are uj^- = l^e 
for K~ and w^o = for K^. 

Calculational details for nucleon phase e^,^^) 
and the K~ phase e~, /i^, K^) have been discussed 
in For the phase {n,p, e~ ,K~ ,K^), with the 
solution of the nucleon and the K~ phase in hand, we can 
calculate energy from Eq. , which keeps decreasing 
as we increase the density. When the condition ujj^o = 
is first achieved, the will occupy a small fraction of 
the total volume and the corresponding charge density, 
qj^a = 0. We can calculate tr, Vq, Rq, kf^ , fc/n, ^/ej 
fc/Ai) PK- and p^o, as in the K~ phase [2J|, with the 
condition w^o = for any chosen baryon density. Af- 
ter obtaining this solution we can calculate the energy 
density and pressure for the K'^ phase. In our earlier cal- 
culations [2J] with higher order couplings, we observed 
that the transition from the nucleon phase to the K~ 
phase is second order in nature. In such transitions the 
mixed phase, where the considered constituents can form 
a cluster, is not favoured. 

For calculations using the E-RMF model with the in- 
clusion of antikaons, we need two distinct sets of coupling 
constants: one being the nucleon-meson coupling con- 
stants and the other, the kaon-meson coupling constants. 
In the former case, the coupling constants are obtained 
by fitting to several properties of finite nuclei [3l| . In this 
work we consider the parameter sets Gl, G2 and FSU2.1. 
The FSU2.1 parameter set [l^ has the same parameters 
as in the FSUGold parameter set [1^ but has one ex- 



tra term in the expression for pressure. The detailed list 
of parameters, for both kaon-meson and nucleon-meson 
couplings (G1,G2 and FSUGold), are given in 



III. RESULTS AND DISCUSSION 

First, we look into the conditions at which the an- 
tikaons start to appear in the symmetric matter, neutron 
matter and NS matter. For this we calculate the antikaon 
energies {ui^- , ^ko) along with difference in chemical po- 
tentials of neutrons and protons (/x„ — Hp), corresponding 
to the considered matter. In Fig. [I] for a given matter, 
the point where ui^- crosses /x„ — jip and the point at 
which ujj^o = represent the onset of K~ and respec- 
tively. These calculations are with the G2 parameter set 
and Uk = -120 MeV. From Fig. [H it is evident that at 
a particular density, lok- increases while w^fo decreases 
with the increase in neutron fraction. This difference is 
purely due to the contribution from the p field [Eq. 
and hence with a large neutron excess the K~ conden- 
sation is less favoured when compared to . However, 
in neutron-rich matter, — /ip is larger and stiffer with 
density. This allows to condense at lower densities. 

Apart from the usual dependence on Uk, the onset 
of condensation of antikaons strongly depends on the pa- 
rameters of the Lagrangian especially for the higher order 
couplings [l^. In case oi K~ , this is due to the strong 
variation in the density dependence of ojk- and /i„ — /^p 
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FIG. 1: Density dependance of the energies of antikaons and 
the difference of the calculated chemical potentials for the 
neutron and proton — Mp) for symmetric matter, neutron 
matter and NS matter. For a given matter, the point where 
ojjf- crosses fin — l-ip and the point at which u^o = represent 
the onset of and respectively. These calculations are 
done with the parameter set G2 with Uk ~ —120 MeV. 



whose interplay determines the onset of condensa- 
tion. The density dependence of w^- is similar to that 
of the EoS and /!„ — /ip varies in a way similar to the 
symmetry energy. So any change in the density depen- 
dence of EoS or that of the symmetry energy will affect 
the onset as well as the effect of K~ condensation. In 
general, K'^ can appear only at densities higher than the 
one at which K'^ condenses [l^. This is due to the fact 
that only with K~ condensation the proton population 
increases and with increasing density the matter becomes 
symmetric at a point where K'^ can start to contribute. 
Thus the onset of K° depends on the onset of which 
in turn depends on the value of Uk and the higher order 
couplings. 

With the G2 parameter set. the densities at which K~ 
and K'^ sets in are (i) for neutron matter : 0.36 and 0.71 
fm-^, (ii) for NS matter: 0.59 and 0.98 fm"^ and (iu) 
for symmetric matter: 1.5 and 1.5 fm~'^, respectively. 
In further discussions, we consider NS matter with Gl, 
G2 and FSU2.1 parameter sets which represent different 
extensions of the RMF model. 

In Figure[51 we present the scalar, vector and isovector 
fields along with the electron chemical potential as a func- 
tion of density for NS matter calculated with Uk = —120 
MeV and —160 MeV. The a and uj fields are attractive for 
antikaons and their pattern is almost similar for all pa- 
rameter sets. The p field is weaker in the case of FSU2.1 
[Figs. [2] (c) and (f )] and it is almost constant (due to the 
additional coupling representing the strength of isoscalar- 
isovector mixing). As discussed earlier, the change in 
energies of antikaons from their symmetric matter value 
depends purely on the p field. Hence the energies of an- 
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FIG. 2: The density dependence of the scalar (qctko'), vector (gujiiVo), and iso-vector {—gpKRo) fields in the NS matter inclusive 
of antikaon phase calculated with Gl, G2 and FSU2.1 parameter sets for Uk ~ —120 MeV (top row) and —160 MeV (bottom 
row). The variation of electron chemical potential also is shown. 
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FIG. 3: The population of hadrons and leptons in NS matter as a function of baryon density calculated using different parameter 
sets (Gl, G2, FSU2.1 shown in columns) with the inclusion of K~ (top row) and both K~ and A'" (bottom row) at Uk ~ —160 
MeV. 



tikaons in the case of FSU2.1 arc not far from their sym- 
metric matter values which will suppress the condensa- 
tion of antikaons, especially the K*^ . Due to this reason, 
with FSU2.1, the condensation happens at a larger 



density and the does not appear at densities relevant 
to NS. In results with the Gl [Figs.[l](a) and (d)], and G2 
[Figs. Hfb) and (e)] parameter sets, as we increase den- 
sity there are two kinks in /ig. The first and second kinks 
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represent the onset of K~ and condensation, respec- 
tively. After the first kink, the p field is enhanced due to 
the increased proton population with the n p + -ft'" 
process. After the second kink, the p field is suppressed 
due to the equal population of protons and neutrons with 
the N N+K^ process which also suppresses the differ- 
ence in and K'^ populations [Eq. ([5])]. The sharpness 
in the kinks increase with the increases in the optical po- 
tential which also leads to the antikaon condensations at 
smaller densities while comparing Figs.[5Ja) and (b) with 
Figs.[2Id) and (e), respectively. 

In Figure |31 the population of different particles in NS 
is plotted against the baryon density, where the calcu- 
lations are done using the GI, G2 and FSU2.1 parame- 
ters with optical potential Uk = —160 MeV. Results are 
given in the presence of only [Figs. [3Ua)-(c)] and with 
both K- and K° [Figs. [3i;d)-(f)]. These results reflect all 
the features that we have discussed with the variation of 
different fields (Fig.[2|). Calculations with Gl, G2 param- 
eter sets in the absence of lead to a situation where 
the proton population exceeds that of neutron population 
[Figs. ISJa) and (b)]. This scenario changes once we in- 
clude the K° in our calculation [Figs.^d) and (e)], which 
suggests that the K'^ starts to contribute at a density 
where the matter becomes symmetric (pp = p„). These 
results predict a symmetric NS core, where the popula- 
tion of both antikaons is almost same. Interestingly at 
higher density, the population of protons (or neutrons) 
and K~ are exactly the same. This means the negative 
charge is solely due to the K~ and hence no leptons are 
present. In case of our results with the FSU2.1 parameter 
set [Figs.[3ljc) and (f)], K~ condenses at a higher density 
and hence A'*^ does not appear, even with Uk = —160 
MeV. The delayed onset of K~ is due to the additional 
coupling as discussed earlier. 

The calculated pressure versus energy density (EoS) 
for different cases is displayed in Fig. H] where we observe 
the regular feature of exotic particles softening the EoS. 
With the onset of K~ the EoS become softer which is 
further softened with the onset of K^. We observe that 
the EoS is strongly influenced by K~ whereas the is 
important only at very high Uk- In the nucleon phase, 
G2 gives a softer EoS compared to Gl and FSU2.1 pa- 
rameter sets. The EoS follows the same pattern after the 
inclusion of K~ and then JsT", with both Uk = —120 and 
— 160 MeV. The sensitivity of EoS in the presence of an- 
tikaons, to the parameter Uk depends on the stiffness of 
symmetry energy. FSU2.1 has a softer symmetry energy 
and hence the corresponding EoS is not very sensitive to 
Uk- 

The role of antikaons in modifying the EoS is very 
well-reflected in the results for the mass-radius relation 
of NS which are presented in Fig. [5] for different param- 
eter sets at different values of Uk- These results are ob- 
tained by solving the well-known Tolman-Oppenheimer- 
Volkoff (TOV) equations [H, [s^ . It is not always true 
that 'earlier is the onset of antikaons, the greater their 
effect on the mass-radius relation" . The different factors 




FIG. 4: The EoS for a nucleon (solid lines), K~ (dashed 
lines) and (dotted lines) phases with (a) Gl, (b) G2, and 
(c) FSU2.1 parameter sets. 



governing the change in maximum mass in the presence 
of A'" are discussed in our previous work [23|. With the 
fact that the presence of K'^ depends mostly on that of 
A'", the influence of K'^ on the mass-radius relation de- 
pends on the influence of AT". Maximum mass with the 
G2 parameter set decreases with the inclusion of K~ at 
Uk = -120 MeV, where K° does not contribute. With 
Uk = —160 MeV, there is a significant reduction of mass 
in the presence of K~ and the inclusion of AT*^ marginally 
decreases the mass further. We observe a similar trend in 
the case of Gl where both the antikaons play a stronger 
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role. With the FSU2.1 parameter set, the sensitivity of 
Uk to the EoS is less and hence the mass and radius 
changes marginally with the inclusion of K~ and the K'^ 
does not play any role. All our results are quite con- 
sistent with the recent observations depicted in Fig. [5] 
It is interesting to note that the recently observed [a| 
pulsar PSR J1614-2230 of mass 1.97 ± O.O4M0 could be 
explained with three different compositions, namely (i) 
with both antikaons in case of Gl, (ii) with K~ only in 
case of FSU2.1 and (iii) without antikaons in case of G2. 
It could be worthwhile to look into the details of these 
compositions. 

Figure |6] represents the population of different parti- 
cles in the NS versus the radial distance from the center 
of the NS, calculated with different parameter sets, when 
the NS has a maximum mass of approximately 2Mq. We 
observe that without antikaons, the NS has a neutron- 
rich core and the presence oi K~ makes it symmetric. 
With the onset of K~ , the population of leptons decrease 
drastically and this change depends on the value of Uk- 
In the case of FSU2.1, we need a large Uk which results 
in a negligible population of leptons in the core of the 
NS. Another interesting feature is the variation of radius 
of a 2Mq NS in the presence of antikaons. To accom- 
modate such exotic particles at a fixed maximum mass 
(M), one should start with an EoS which is stiffer and 
yields a larger maximum mass without exotic particles. 
The presence of an exotic core reduces the maximum 
mass whereas the corresponding radius (R) may cither 
increase or decrease. The antikaons can increase or de- 
crease the central baryon density (pc) [11] and hence R 
(i? oc l/pc), however only marginally. To accommodate 
more antikaons we need to start with a stiffer EoS yield- 
ing larger M and R, without antikaons. Thus for a given 
maximum mass, the corresponding radius is more if we 
have more antikaons (exotic particles). However, we may 
not observe stars extremely close to the maximum mass 
and hence the corresponding radius could vary. Hence 
more work is needed to understand the sensitivity of the 
radius of massive NS to the presence of exotic particles. 
In this work, we have ignored the presence of hyperons 
that can affect the role of antikaons and the NS proper- 
ties. The onset of some of hyperons (A and E") can be 
in the same density range corresponding to the onset of 
antikaons. Hence the presence of hyperons can not only 
soften the EoS further but could cause a strong inter- 



play between the hyperons and antikaons, affecting each 
other's role on the EoS. It will be interesting to study 
all these effects within the extended RMF models with 
higher order couplings. Work is in progress in this direc- 
tion. 



IV. CONCLUSIONS 

With the inclusion of both antikaons {K~ and 
in extended relativistic mean-field models (with param- 
eter sets Gl, G2 and FSU2.1), we observe that the on- 
set of condensation of antikaons strongly depends on the 
kaon optical potential {Uk) and the parameters of the 
Lagrangian, especially the higher order couplings. This 
is similar to the conclusion in our earlier work [2J] done 
only with the inclusion of , where we attributed the 
onset as well as the effect of K~ condensation to the 
change in the density dependence of the EoS or that of 
symmetry energy. As can appear only at densities 
higher than the one at which K~ condenses, the onset of 
K° depends on the onset of K~. With Gl and G2, we ob- 
serve that the EoS is strongly influenced by K~ whereas 
the K° is important only at higher Uk (> -140 MeV). In 
case of the FSU2.1 parameter set, the additional higher 
order coupling softens the symmetry energy and hence 
the K~ condensation happens at a larger density and 
the K'^ does not appear at densities relevant to the NS. 
These effects are well-reflected in the mass-radius rela- 
tion and the composition of NS. The onset of K'^ leads to 
symmetric and lepton-deficient matter at the core of the 
NS which would be proton-rich if we ignore K° while the 
K~ dominates. We also observe that a 2Mq NS can be 
explained in three ways with: (i) a stiffer EoS with both 
antikaons, (ii) a relatively soft EoS with K~ and (iii) a 
softer EoS without antikaons. In the case of 2Mq being 
the maximum mass of an NS, we observe that greater 
concentration of antikaons leads to an increase in the ra- 
dius (R-i) of such stars. Without antikaons (G2) we get 
i?2 = 11.03 km, with K- (F_SU2.1) we get i?2 = 11-42 
km and with both K' and K° (Gl) we get i?2 = 12.55 
km. It would be interesting to study whether a precise 
information about the radius of massive NS could reveal 
the presence of exotic cores. 
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FIG. 5: The mass-radius relation from extended RMF models. 
Different curves represent the calculations using different kaon 
optical potentials (Uk) and with different parameter sets. For 
each parameter set, solid black line represents the pure nu- 
cleon phase, lines with different patterns and colors represent 
the phase with K~ and the corresponding small dotted lines 
represent the phase with both antikaons (K'jK"). The dif- 
ferent patterns and colors represent the strength of the kaon 
optical potential Uk {\Uk\ quantifies the influence of kaons) 
as specified in the inset. The solid circles represent the maxi- 
mum mass in every case. Mass is given in units of solar mass 
Mq. Solid squares {vph = R) and open triangles (r^h ^ R) 
represent the observational constraints 0|, where rph is the 
photospheric radius. The shaded region corresponds to the 
recent observation of a f .97 ± O.O4Af0 star 
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FIG. 6: The composition of an NS with a maximum mass of 
2Mq, calculated using G2, FSU2.1 and Gl parameter sets. 
The number density of various particles is plotted against the 
distance from the center of the NS. The values of Uk are 
adjusted to have the desired maximum mass. 



